Metal-insulator transition in 2D random fermion systems of chiral symmetry classes 
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Field-theoretical approach to Anderson localization in 2D disordered fermionic systems of chiral 
symmetry classes (BDI, AIII, CII) is developed. Important representatives of these symmetry 
classes are random hopping models on bipartite lattices at the band center. As was found by Gade 
and Wegner two decades ago within the sigma-model formalism, quantum interference effects in 
these classes are absent to all orders of perturbation theory. We demonstrate that the quantum 
localization effects emerge when the theory is treated non-perturbatively. Specifically, they are 
controlled by topological vortex-like excitations of the sigma models. We derive renormalization- 
group equations including these non-perturbative contributions. Analyzing them, we find that the 
2D disordered systems of chiral classes undergo a metal-insulator transition driven by topologically 
induced Anderson localization. We also show that the Wess-Zumino and Z2 theta terms on surfaces 
of 3D topological insulators (in classes AIII and CII, respectively) overpower the vortex-induced 
localization. 

PACS numbers: 73.23.-b, 73.20.Fz, 73.43.Nq 



I. INTRODUCTION 

Disorder-induced metal-insulator transitions represent 
a fundamental concept of modern condensed matter 
physics. The principal mechanism of this phenomenon 
is Anderson localization 1 - due to quantum interference 
between different trajectories of electron propagation in 
a disordered medium. As was realized in seventies by 
Thouless and Wegner^ 3 -, the problem of Anderson lo- 
calization is closely connected to the scaling theory of 
critical phenomena. Several years later, Abrahams and 
coauthors formulated the scaling theory of localization 4 
describing the problem in terms of a flow of the dimen- 
sionless conductance g with the system size L. This 
phenomenological theory was put on a solid basis by 
Wegner— who discovered the field-theoretical description 
of Anderson localization in terms of a nonlinear sigma 
model. The latter is an effective field theory of the 
problem, capturing all essential properties of the sys- 
tem including its symmetry and topology. The crucial 
importance of the topological aspect was first demon- 
strated by Pruisken in the framework of the quantum 
Hall effect^. The field-theoretical (sigma model) formu- 
lation of the problem underpins the notion of universal- 
ity of critical behavior at metal-insulator transitions, i.e., 
its independence of microscopic details of the electron 
motion. Over the years the sigma model permitted the 
great progress in understanding localization phenomena 
including, in particular, critical behavior ^Si 8 . energy level 
and wave function statistics;^ symmetry classification^ 
interaction (most prominently 2D metal-insulator transi- 
tion) and non-equilibrium effects* 1 - 1 - - — theory of uncon- 
ventional disordered superconductors,— classification of 
topological insulators and superconductors ^ 

Recent years have witnessed increasing interest and 



significant progress in the field of the Anderson local- 
ization. Development of the complete symmetry classifi- 
cation of disordered systems 1 ^ is one of most important 
advances. It has been realized that underlying symme- 
tries and topologies induce a remarkable variety of the 
types of infrared (long-distance, low-temperature) behav- 
ior, including, in particular, critical phases and quantum 
phase transitions between metallic and insulating states 5 . 
Experimental discoveries of graphene 1 ^ and of topologi- 
cal insulators 1 ^— have given an additional boost to the 
field. Disordered Dirac fermions in graphene in the ab- 
sence of valley mixing and surface states of disordered 
topological insulators and superconductors are character- 
ized by sigma models with Wess-Zumino or theta terms, 
ensuring a topological protection from localization ^ 16 ' 21 . 

The full symmetry classification of disordered fermion 
systems&iS includes three families of symmetry classes: 
conventional (Wigner-Dyson), chiral, and superconduct- 
ing (Bogoliubov-de Gennes). In this paper we will con- 
sider localization properties of two-dimensional (2D) chi- 
ral models. Characteristic feature of these systems is 
the chiral symmetry of the underlying Hamiltonian H 

('il< where C 2 = 1. It is convenient to use the 
block structure such that C is represented by the third 
Pauli matrix r z ; then the Hamiltonian acquires a block 
off-diagonal form, 



H = 



h 
h) 



(1) 



One-dimensional disordered systems with chiral sym- 
metry were studied in Refs. [22|. An effective field- 
theoretical description - non-linear sigma model - for 
chiral Hamiltonians with randomness was developed in 
Refs. I23U241 Spectral properties of chiral random matrix 
ensembles were discussed in the context of mesoscopic 
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transport and quantum chromodynamics in Refs. [25ll26l . 

A standard realization of a chiral system is provided 
by a bipartite lattice model (i.e., containing two sublat- 
tices) with random hopping between the sublattices and 
no diagonal disorder. The Hamiltonian of such a model 
at the band center is manifestly block off-diagonal in the 
sublattice space. In analogy with Wigner-Dyson classes, 
there are three chiral classes: chiral unitary (AIII), chi- 
ral orthogonal (BDI), and chiral symplectic (CII) classes. 
The class AIII describes chiral systems with broken time- 
reversal symmetry, which may be due to magnetic flux 
(random or homogeneous) through the lattice plaque- 
ttes. The classes BDI and CII are characterized, in addi- 
tion to the chiral symmetry, by time-reversal invariance, 
H = KH K^ 1 , where the unitary matrix K satisfies 
K 2 = ±1. For the chiral orthogonal class BDI one has 
K 2 = 1; in the standard representation K = 1. For 
the chiral symplectic class CII, which corresponds to the 
case of broken spin rotation symmetry (usually due to 
spin-orbit interaction), K 2 = —1; in the standard repre- 
sentation K = ia y , where a y is the second Pauli matrix 
in the spin space. Thus, in the standard form, random 
Hamiltonians of the chiral classes have block off-diagonal 
structure and consist of complex, real, and real quater- 
nion entries for the classes AIII, BDI, and CII, respec- 
tively. 

Realizations of 2D chiral classes arise naturally when 
one models the disordered graphene at the Dirac point^l. 
Specifically, the chiral symmetry occurs if the dominant 
disorder is of random-hopping character. An important 
example is the random magnetic flux, either genuine or 
effective (due to lattice corrugations), in which case the 
system belongs to the chiral unitary class AIII. Another 
possible realization is provided by graphene with vacan- 
cies. A vacancy can be modeled by cutting all lattice 
bonds adjacent to the vacated site, so that vacancies rep- 
resent a special type of bond disorder. This model is in 
the chiral orthogonal symmetry class BDI. If an external 
magnetic field is applied to graphene with vacancies, the 
time-reversal symmetry is broken and the system falls 
into the class AIII. Taking into account the spin-orbit 
interaction (in the absence of magnetic field) brings the 
system into the chiral symplectic symmetry class CII. 

Clearly, the above symmetry analysis is not specific for 
graphene (hexagonal lattice) but is equally applicable to 
a random-hopping model on any other bipartite lattice 
(e.g. the square lattice) at the band center. The analysis 
of Anderson localization in the chiral classes performed in 
the present work is equally applicable to all such models. 

A remarkable property of the three chiral symmetry 
classes in 2D is the exact absence of weak localization 
corrections to all orders in the perturbation theory^. At 
the same time, when the disorder is strong enough, the 
transition to insulating state is inevitable. Indeed, con- 
sider a tight-binding Hamiltonian on a 2D bipartite (e.g., 
square or honeycomb) lattice. Let us introduce disorder 
by randomly cutting the lattice bonds. (As discussed 
above, this model belongs to the chiral orthogonal class 



BDI.) Classically, this is the standard percolation model 
with bond disorder. When the concentration of removed 
bonds reaches the critical value (1/2 for square lattice 
and 1 — 2sin(7r/18) (=> 0.65 for honeycomb lattice), the 
system becomes disconnected and undergoes the classi- 
cal percolation transition to the insulating state. Obvi- 
ously, the conductivity of such a system (defined in the 
limit of the infinite system size) is zero. It is well known 
that, in Wigner-Dyson symmetry classes, such a classical 
percolation transition is preceded by the Anderson metal- 
insulator transition. Specifically, due to quantum inter- 
ference, the conductivity of a quantum system becomes 
zero at the point where classically it would still conduct 
(the remaining bonds are still percolating). Does this 
happen also for models of chiral classes? The exact ab- 
sence of quantum-interference effects to all orders^ sug- 
gests that the answer is negative. However, numerical 
I II /-J 

simulations in Refs. 28.29 do show a quantum localization 
transition in chiral classes. These papers explored numer- 
ically models with continuous disorder which are always 
conducting classically, and yet found a metal-insulator 
transition. Thus, there appears to be a contradiction 
between the absence of localization in the sigma-model 
analysis of Ref. [24] and the numerical results. This con- 
tradiction is resolved in the present paper. 

We will show below that, while the analysis of Ref. 
l24l is fully correct to all orders of perturbation theory, 
the quantum interference effects do emerge in the sigma- 
model field theory when one considers non-perturbative 
effects related to topologically non-trivial, vortex-like ex- 
citations. The field-theoretical analysis of the quantum 
localization in chiral classes requiring a systematic study 
of such non-perturbative effects constitutes the subject 
of the paper. 

Non-perturbative effects caused by vortices were con- 
sidered in the works of Berezinskii^S and Kosterlitz and 
Thoulcss 31 in the context of planar 2D ferromagnet (XY 
model) and 2D superfluid flow. These systems are char- 
acterized by a U(l) field (direction of magnetization and 
phase of the condensate wa)ve function, respectively) and 
exhibit the temperature-driven Berczinskii-Kosterlitz- 
Thoulcss (BKT) phase transition between quasi-ordered 
(ferromagnetic or superfluid) and disordered (paramag- 
netic or normal fluid) phases caused by vortex unbind- 
ing. Renormalization group (RG) analysis of this phase 
transition was first developed by Kosterlitz^. Alterna- 
tive descriptions of the same systems are provided by the 
Coulomb gas (see Ref. [H for a review) and sine-Gordon 
models. RG equations for the sine-Gordon model, equiv- 
alent to Kosterlitz RG, were derived in Ref. [3J in the 
one-loop order and later in Ref. HH with the two-loop 
accuracy. We will show that, in analogy with the BKT 
physics, vortices drive the localization transition in chiral 
classes. However, the emerging theory is more compli- 
cated and the character of transition is essentially differ- 
ent, since the symmetry group of the localization problem 
is different from U(l). 

The layout of the paper is as follows. In Sec. [Hi we 
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present the general structure of the sigma model for chi- 
ral symmetry classes, outline the proof of the absence of 
perturbative weak localization corrections, and identify 
the mechanism of non-pcrturbativc effects. The general 
background-field formalism for the renormalization of the 
sigma model is developed in Sec. IIIII and then applied 
to the calculation of perturbative corrections in Sec. IIVI 
In Sec. El we construct the vortex excitations and ex- 
plore their effect on the conductivity. For this purpose, 
we derive a set of RG equations on the conductivity <7, 
Gade-term coupling c (controlling the behavior of the 
density of states), and the vortex fugacity y. Its analysis 
reveals a metal-insulator transition. Therefore, we show 
that the chiral classes do exhibit quantum interference 
effects - the weak localization and the Anderson tran- 
sition - driven by topological excitations in the sigma 
model. In Sec. IVII we develop an alternative approach to 
non-perturbative renormalization in class AIII. This ap- 
proach is based on a mapping of the original sigma-modcl 
with vortices onto a dual theory which has a character of 
generalized sine-Gordon theory. Section IVIII is devoted 
to the interplay of this topological localization with fur- 
ther possible topological properties of 2D chiral-class sys- 
tems: the Wess-Zumino term in class AIII and the Z2 
theta term in class CII. These terms in the sigma models 
emerge when the 2D systems are realized on surfaces of 
3D topological insulators of the corresponding symme- 
try classes. We show that the vortex-induced topological 
localization becomes inefficient in the presence of such 
terms. The latter thus ensure the topological protection 
of the system from localization. Our results are summa- 
rized in Sec. I Villi Some technical details are presented 
in two Appendices. 



II. PRELIMINARIES 

Very generally, localization properties of a disordered 
system are described by the non-linear sigma model field 
theory. We will use the fermionic replica version of 
the sigma model. (Alternatively, the same analysis can 
be performed within the supersymmetry approach; we 
choose the replica formalism here since in this case the 
presentation turns out to be physically more transpar- 
ent.) The sigma-model target space is a compact sym- 
metric space, see Ref. [H for a review. For the three chiral 
classes, the field Q(r) is a unitary matrix taking values 
in the symmetric spaces listed in Table |TJ The size of the 
Q matrix, N, is determined by the number of replicas; in 
the end one should take the replica limit N — > 0. In each 
of the classes CII and BDI, the Q matrix obeys a certain 
linear constraint. These constraints can be represented 
as Q = Q, where the "bar" operation is defined in the 
following way: 

Q-|° T ' CI1 ' (2) 
^-\CQ T C, BDI. [) 



Class 


Manifold 


Q matrix 


D 


s 


AIII 


XJ(N) 


No restriction 


N 2 


1 


CII 


U(N)/0(N) 


Q = Q T 


N(N + l)/2 


2 


BDI 


U(JV)/Sp(JV) 


Q = CQ T C 


N(N - l)/2 


2 



TABLE I: List of chiral symmetry classes including corre- 
sponding compact sigma-model manifolds, symmetry con- 
straints on the Q matrix, and the dimensionalities D of the 
target spaces. Parameter s in the last column is introduced 
to normalize the conductivity in Eq. <(3j . 



In the class BDI, C is a skew-symmetric matrix such that 
C 2 = 1. In the standard parametrization of the model, 
the number of replicas N is even and includes the spin 
degree of freedom while the matrix C is the a y Pauli 
matrix acting in the spin subspace. 

The distinguishing feature of the three manifolds of 
chiral-class sigma models is that the corresponding tan- 
gent spaces possess a unit generator commuting with all 
other (traceless) generators. This generator corresponds 
to the U(l) part of the corresponding symmetric space. 
The presence of this special generator allows an addi- 
tional Gade term in the sigma model action 2 ^: 

S[Q] = - J ^ [a TriQ- 1 VQ) 2 + c(Tr Q~ l VQf] . 

The parameter a is the dimensionless conductivity of the 
sample (measured in units e 2 /irh)&nd c multiplies the 
Gade term. The bare value of c is of order 1. 

Let us decompose Q according to Q = e % ^U with 
det U = 1. The overall phase is exactly the U(l) degree 
of freedom characteristic for sigma models of the chiral 
symmetry classes. The action Q takes the form 

S[Q]= [ ^1 \N(a + Ac)(V0) 2 + a TrfVC/- 1 VC7)1 . 

J 87TS 

(4) 

We see that variables cj> an d U decouple and that the 
action is quadratic in <f>. This means that the prefac- 
tor a + Nc of the (V0) 2 term is not renormalized. In 
the replica limit N — > 0, the absence of localization cor- 
rections to a follows, which is the result of Gade and 
Wegner—. 

The above argument for the absence of the renormal- 
ization of conductivity has, however, a caveat. Specifi- 
cally, the action for the variable <fi is not strictly Gaussian 
in view of the compact U(l) nature of this degree of free- 
dom. The group U(l) is not simply connected allowing 
for topological excitations - vortices. Fluctuations of the 
matrix Q involving vortices yield corrections to a + Nc 
via the BK T 30 ' 31 mechanism. These corrections break 
the replica symmetry and thus invalidate the Gade and 
Wegner argument. Our goal is to analyze the effect of 
vortices on the renormalization of the chiral-class sigma 
models. We will start with the general background field 
formulation of the renormalization. Then we will con- 
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sider both perturbative and non-perturbing (due to vor- 
tices) contributions to the RG equations. 



III. BACKGROUND FIELD FORMALISM 

Very generally, renormalization of the sigma model can 
be implemented in the framework of the background field 
formalis m 36 ' 37 . Following this approach, we first separate 
slow and fast parts in the matrix Q. The former is re- 
ferred to as the background field. Then the partition 
function is integrated over the fast part of Q yielding the 
effective action for the background field. This effective 
action has the same sigma-model form with the renor- 
malized parameters a and c. 

In order to renormalize the action ([3]), we parametrize 
Q — VQV, where Q is the fast field of the same symmetry 
as Q and unitary matrices V and V represent the slow 
background field. In the symmetry class AIII, matrices 
V and V are independent while for the other two chiral 
classes they are related by Eq. @ . After integrating out 
Q, the effective action involves only the gauge invariant 
combination of slow fields Q' — VV. It is convenient to 
introduce the following notations for the gradients of the 
slow field: 



between the slow and fast modes read 



A = WV~ 



A = V- X VV. 



(5) 



Matrices A and A are related by Eq. @ in classes CII 
and BDI. 

We start the renormalization group analysis with the 
action Sq[Q] given by Eq. with bare parameters <To 
and Co- Upon substitution Q — VQV, the action decom- 
poses into four terms: 



So [Q] = So [Q'\ + S L [Q] + s£> [V, V, Q] + S^> [V, V, Q] . 

(6) 

The first term, So[Q'], is the bare action for the slow 
background field Q' — VV. It can be represented as 

^£ {a Tr (A + A) 2 + c [Tr (A + A)] 2 } . 

(7) 

The action for the fast field Sl[Q) is given by the same 
Eq. ([3]) with bare parameters <Jo and cq but includes also 
the mass term, 



:( 2 )r 



S L [Q] =S [Q] 



(TO 



8irsL 2 



d 2 rT r (2-Q-Q- 1 ). (8) 



The mass ~ 1/L is introduced to ensure that the ma- 
trix Q contains only fast degrees of freedom. At the 
same time, we will assume that the slow fields V and 
V change very little on the distances of the order of L. 
It is important that the mass term is symmetric on the 
manifold of Q; this will guarantee that the renormalized 
action involves only the gauge-invariant combination VV 
and retains the full symmetry of the sigma model. Fi- 
nally, the terms in the action representing the interaction 



S 4 



(1) _ 



a r 



go Tr VQQ^A + Q- 1 VQA 



c Tr Q~ l \7Q Tr (A + A) 



S' 1 



(2) 



AA - Q _1 AQA 



(9) 
(10) 



Integrating e~ s ^ over fast modes Q, taking the log- 
arithm of the result, performing the gradient expansion, 
and retaining the terms with two spatial derivatives (i.e. 
up to second order in A and A), we obtain the effective 
action for the slow field: 



S[Q'} = S [Q'] 



-42) 



(1)- 

int _ 



(ii) 



Here (...) denotes the averaging over Q with the weight 
e -s L [Q}_ Note that the average (S^) is linear in the 
gradients of the slow variables and is hence zero. The 
renormalizability of the sigma model implies that the ef- 
fective action for the slow field Q' has the form of Eq. ([7]) 
with the bare parameters <7o and cq replaced with their 
renormalized values a and c. Using this fact, we will as- 
sume three particular forms of the background field and 
thus establish the renormalized parameters in terms of Q 
integrals for all three symmetry classes. 

The D-dimensional sigma-model target space is gen- 
erated by D — 1 traceless Hermitian generators and one 
generator proportional to the unit matrix. In the classes 
CII and BDI these generators are subject to an additional 
symmetry constraint T — T . We normalize the genera- 
tors as Tr T a T b = S ab . The traceless generators of the 
sigma-model manifold obey the following Fierz identities: 



> 3«Iw 



1 



1 

N 



• Sad, 



ijOkl, 



Sik$ji) 



' 1 



- {8n6jk + CikCji) ~ — SijS. 



AIII, 



CII, 



BDI. 



(12) 

Assume the background field of the form A = A = iJT 
with some constant vector J and T being one of the D — 1 
traceless generators of the sigma-model manifold. Sub- 
stituting this background field into Eq. (fTTj) . averaging 
over directions of J, and comparing the prefactors of J 2 
in the left- and right-hand side of the equation, we find 
the renormalized conductivity in the form 



°o-^Jd*rTr([Q,T][Q-\T] 



32irsA 



d 2 rTr{VQ,Q _1 }r 



1 2 ' 



(13) 



Next, we assume a pure gauge background field such 
that A = — A = iJt. The matrix t is now one of gen- 
erators of the stabilizer group, U(A r ), O(iV), and Sp(7V), 
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for the classes AIII, CII, and BDI, respectively. These 
generators obey the second Fierz identity 



'SuSjk, 
2 (SuSjk — CikCji) 



AIII, 

CII, 

BDI. 



(14) 



Such a gauge background field corresponds to the con- 
stant Q' = 1 and hence the left-hand side of Eq. (fTTj) 
vanishes. Averaging over directions of J, we obtain the 
following Ward identity 



° = ||/ d 2 rTr([Q,t] [Q'\t] 



(15) 



We sum up identities (TTBI with all generators T and 
add identities (IT5|) with all possible matrices t. With the 
help of the Fierz identities (|12[) and (fT4"]) and using the 
properties Q = Q in the classes CII and BDI, we obtain 
the renormalized conductivity 



By 



(16) 



Here we have introduced the averages 



L>\ = I ^^(^(Q-^QUQ-^Q)^), (17) 



B 2 = 



8ns 
d(r - r') 
8ns 



(18) 



Finally, consider the background field of the form A = 
A = iJl/y/N generated by the unit matrix. This field 
configuration brings into play both the kinetic and the 
Gade terms of the action. After averaging over directions 
of J, we obtain the identity 



Nc- 



CO 



AT , (c + Acq) 2 

Nc H — B 2 . (19) 



With a from Eq. (|TB)) . this yields the following expression 
for the renormalized coupling c: 



Co 



N(D-l) 



D 



Equations (116)) and (|20|) represent renormalization of 
a and c in the most general form. The averages B\^ 2 
implicitly depend on the infrared scale L set by the mass 
term in the fast field action ©. The RG flow equations 
are obtained in the standard way by taking the derivative 
of Eqs. (fTl))) and (|20l) with respect to InL and replacing 
do a- n d Co in the right-hand side of these equations with 
the running, L-dependent values a and c of the couplings. 



IV. PERTURB ATIVE RENORMALIZATION 

In the previous section we have developed the back- 
ground field formalism and reduced the renormalization 
of the sigma models of chiral classes to the calculation of 
the correlators (|17p and (fT5|) . The matrix Q belongs to 
the corresponding sigma-model manifold and averaging is 
performed with the statistical weight determined by the 
action ([8]). We will now perform this averaging using a 
saddle point method. The saddle point approximation is 
justified in the limit oq 3> 1. We start with the spatially 
uniform saddle point Q = 1. Note that the mass term in 
the action ([S]) is minimized by this, rather than any other, 
constant value of Q. Expansion in the vicinity of the 
spatially uniform saddle point yields perturbative contri- 
butions to the renormalized couplings. As we are going 
to explain, such perturbative correction to a vanishes, in 
the replica limit N — > 0, in all orders of the perturbation 
theory for all three chiral classes. We will then proceed 
by including a non-perturbative contribution from saddle 
configurations containing a vortex-antivortex dipole. We 
will show that this yields a non-zero result. 

We parametrize the fast field Q = 1 + iW - W 2 /2 + 
0(W 3 ) by a Hermitian matrix W subjected to the lin- 
ear constraint W = W in classes CII and BDI. Let us 
decompose W in generators of the corresponding sym- 
metric space: W = wqI/vN + ^2 a w a T a . Expanding 
Eq. (JSJ to the second order in W, we obtain the Gaus- 
sian action 



S L [W] 



d 2 r 



(a + Ac )(Vu>o) 2 +<J L 2 wl 



C() 



+ L- V 



(21) 



This action yields the following propagators of the com- 
ponents of W: 



(wo(q)wo(-q)) = 



4ns 



(22) 



(era + Nc )q 2 + co^ 2 ' 
4:7rs 

(w a (q)w a (-q)) = - (23) 
er (<r + L J 

Within the perturbative calculation, the average B 2 
is identically zero since it involves only the U(l) com- 
ponent of Q; specifically, TrQ _1 VQ = iy/NVwo. Con- 
sequently, £?2 vanishes, B 2 oc q 2 (wo(q)WQ(— q)) 9 _>o = 0, 
that implies that cr+Nc is not renormalized, see Eq. (|19l) . 
This cancellation of B 2 is the essence of the Gade-Wegner 
argument for the absence of weak localization correc- 
tions to conductivity. We will see later that the saddle 
points involving vortices yield a nonvanishing contribu- 
tion to B 2 . 

Contrary to B 2 , the average B\ is non-zero already 
on the perturbative level. We will calculate it within 
the one-loop approximation. Rewriting Eq. (|17p in the 
momentum representation, expanding in the components 
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of W, and applying the Wick theorem, we find 
q 2 d 2 q 



Bi = 



32tt 3 s 



{W l] (q)W kl (-q)){W jk (q)W H (-c l )) 



(^(q)W- jfc (-q))(^(q)^(-q)) . (24) 



We have dropped the term ~ q 2 (W 2 )\ q ^o- The absence 
of such a term is justified by the finite mass in the action 
(|2"Tj) and hence a finite limit of the W propagator at q —> 
0. Using the decomposition of W in generators T and 
the propagators (|22]) and (|23|) , we obtain B\ in the form 



Bi = 



q 3 dq 



2a 2 (q 2 +L 



-2\2 



^Tr [T a ,T b Y 



(25) 



a . b 



Note that the propagator (wqWq) is canceled in the ex- 
pression for B\. The momentum integral is logarithmi- 
cally divergent and we have introduced the ultraviolet 
cut-off A to regularize this divergence. 

Using the appropriate Fierz identity, we calculate the 
trace of products of four generators and obtain the result 



(26) 



Substituting this result and B2 = into Eqs. ([16]) and 
(f20|) . we obtain logarithmic corrections to a and c that 
can be recast in the form of renormalization group equa- 
tions 



da 
d\nL 

dc 

d\nL 



-N + NO(l/a), 
l + 0(l/a). 



(27) 
(28) 



These equations describe the real space scaling with the 
running infrared cut-off length L. The higher loop con- 
tributions provide 0(1/ a) terms in the perturbative RG 
equations. In the replica limit, renormalization of a van- 
ishes, in full agreement with the argument by Gade and 
Wegner— . 

A remarkable feature of the perturbative renormaliza- 
tion in class AIII is that the right-hand side of Eq. (|2"8)) is 
exactly 1 in the limit N — > 0. Higher loop contributions 
to the beta function for c are proportional to N and van- 
ish in the replica limit . 38 i 39 This is not the case in the 
other two chiral symmetry classes. The exact perturba- 
tive RG equations in class AIII can be applied even when 
a is not large. We will use this fact below for constructing 
the non-perturbative RG flow diagram. 



V. NON-PERTURBATIVE 
RENORMALIZATION 

In order to find a correction to conductivity that does 
not vanish in the limit N — > 0, we have to consider the 
saddle-point configurations other than Q = 1. These 



configurations will include vortices. Vortices are singular 
points of the Q matrix such that the overall phase </> of 
Q rotates by 2n along any path going around the vortex. 

The sigma-model description with small gradients of 
Q breaks down at distances of the order of the mean 
free path from the vortex center. This region, referred 
to as the vortex core, should be excluded from the sigma 
model action The contribution of the vortex core to 
the overall action, S v , is not universal and depends on 
details of the ballistic electron dynamics inside the core. 
Generically, S v oc cto- This introduces a new coupling 
constant in our theory, referred to as fugacity, yo = e~ Sv , 
that is the statistical weight associated with the vortex 
core. 



A. Vortex configurations 

Let us construct explicitly the vortex configuration in 
three chiral classes. The minimal model of classes AIII 
and CII with a single replica reduces to the abelian U(l) 
model with Q = e 1 ^ . The same situation occurs in class 
BDI at the minimum number of replicas N = 2. In the 
latter case the matrix Q retains the 2x2 spin structure: 
Q = e l ^\. In all three models, the action (|5]) has the 
form (we neglect the mass term for simplicity) 



S[Q] 



K 
2^ 



d 2 r(V(t)f 



where we have introduced the stiffness parameter 



K 



(*■ 

IK 



c)/4, 
c)/8, 
2c)/4, 



AIII, 

CII, 

BDI. 



(29) 



(30) 



Assume a set of vortices with positions and charges 
(vorticities) rii = ±1. Such a configuration creates the 
gradient of </> given by 



V <; 



eq&y^»iV&hi|r 



(31) 



Indeed, this field <f> satisfies the equation of motion V 2 </> = 
with additional condition of 2irni winding of (f> around 
every vortex. The mass term in the action which we 
have temporarily neglected in Eq. (|2"9")l , requires <p = 
at infinity. This is only possible if the system of vortices 
is neutral: J^- n; = 0. Thus the simplest allowed con- 
figuration of vortices is a vortex-antivortex pair, i.e., a 
dipole. Configurations with more than two vortices or 
vortices with higher winding numbers have lower statis- 
tical weight and can be safely neglected. 

Let us consider an individual dipole with vortex and 
antivortex at positions rj. 2< The separation between 
these points, m = ri — r2, will be called the dipole mo- 
ment. Substituting Eq. (JSJ) into Eq. (gU), we find the 
action for the dipole: 



2S V +2jnn(A|m|). 



(32) 
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Class 


Null-space 




Vp 


AIII 
CII 
BDI 


cpN -l 
R piV-l 


7T 

~N/2- 
7T 


-Vr(iv) 

Vr(iv/2) 

~ 2 /r(AO 



TABLE II: List of spaces spanned by equivalent dipole con- 
figurations. Volumes of these spaces are denoted by V p . 



Here A is the ultraviolet cut-off of the order of inverse ra- 
dius of the vortex core. The action 2S V associated with 
two vortex cores is explicitly added. We see that the 
action of the dipole grows logarithmically with increas- 
ing m. This is true provided the dipole is smaller than 
the infrared length L. For larger dipoles, the mass term 
will modify the dipole saddle point, and the action will 
grow linearly with m. Thus we conclude that the infrared 
cut-off L in the Q theory effectively limits the maximum 
allowed size of a dipole. 

Let us now consider a general situation with arbitrary 
number of replicas. In class AIII, we place the dipole in 
the first replica and represent the matrix Q as 



= l + |p)(e**-l)<p| 



(33) 



Here \p) = {1,0,0,...} is the iV-component vector in 
the replica space. We can generate other dipole configu- 
rations with the same action (|32[) by moving and rotating 
the dipole in the real space, which changes just the func- 
tion <f>, and also by rotations Q m VQ m V in replica 
space with spatially constant unitary matrices V and V. 
The mass term in Eq. ([8]) requires Q m — 1 at infinity 
and thus restricts the matrices by VV = 1. Such unitary 
rotations lead to the same form of Q m given by Eq. ([55)1 
with an arbitrary iV-dimensional complex unit vector \p) . 
Thus equivalent dipole configurations in replica space 
span the manifold CP^ 1 = U(A)/U(A - 1) X U(l) = 

In class CII, additional constraint V = V T applies, see 
Eq. ([2]). As a result, allowed rotations generate an ar- 
bitrary unit real vector \p) with N components. The 
corresponding space is RP^ -1 = 0(N)/0{N — 1) x 
0(1) = S N ~ 1 /S . Finally, in class BDI, the constraint 
@ leads to the unit vector \p) composed of N/2 real 
quaternions. The corresponding manifold is HP^ 2-1 = 
Sp(A)/S P (A - 2) x Sp(2) = S^^/S 3 . The spaces of 
equivalent dipole configurations are listed in Table [TT] to- 
gether with their volumes V p . In all three chiral symme- 
try classes V p = 0(1/ N) in the replica limit N — >• 0. 



around the dipole saddle point. The result has the form 

Z' 



B^ = A* I al'Rd 



(34) 



Here S m is the action (j32|) associated with the dipole, 
■^dip an d Z-q are the partition functions of quadratic fluc- 
tuations near the dipole configuration and near the trivial 
uniform minimum of the action, respectively. Zero modes 
related to rotations of the dipole in replica space are ex- 
cluded from the determinant Z' dip ; they yield the volume 
V p of the space of equivalent dipoles, see Table [TTJ Zero- 
modes related to the position R of the center of mass of 
the dipole and nearly-zero modes related to the dipole 
size m are excluded from Z' dip as well and are taken into 
account in Eq. (|34[) via integration over the correspond- 
ing collective coordinates. The factor A 4 arises when the 
partition sum over positions of vortex and antivortex, 
with the core size ~ 1/A, is replaced by the integral over 
R and m. (Appearance of this factor is clear on dimen- 
sional reasons.) Finally, the pre-exponential factors 
are given by expressions inside angular brackets in Eqs. 
ffTTJl and CE5]) evaluated on Q m from Eq. ([33]). 

The two vortices forming the dipole are located at R± 
m/2. From Eqs. ([311) and (ESI), we obtain 



Qm ^ uQr, 



-iV^(r) \p)(p\ 



In 



r R m/2 



r- R + m/2 



\p){p\. (35) 



In order to calculate , we substitute this expression 
into Eqs. (JT7J and flg]l. Afterwards, Eq. ([3J]) is used to 



find B^ 1 ^. The projective property of the matrix 
and the trace structure of Eqs. (fTT|) and ([T8|) imply the 
relation 



B. 



dip 



uB 



d i ] :> 



(36) 



where we have introduce the notation u = Tv\p){p\. In 
the classes AIII and CII, u = 1, while in the class BDI 
we have u = 2. The quantity u can be interpreted as the 
minimal value of the number of replicas N allowed by the 
symmetry. 

In view of the identity ([36]). it suffices to calculate B^ 9 . 
When the dipole position is fixed, the integrand in Eq. 
dl~7l) depends explicitly on both r and r' rather than on 
their difference only. The uniformity is restored after 
integration over R in Eq. (|34l) . Thus we will carry out 
the integration over the dipole positions first. 

Equations (fT7[) . (|34)) . and (|35|) combine into 



B. Derivation of RG equations 

The contribution of the dipole to the correlation func- 
tions Bi t 2 of the fast field Q is found by substituting 
Q = Qm from Eq. into Eqs. ([TTJ) and (JJ5J) and cal- 
culating the Gaussian integral of quadratic fluctuations 



B'l 



dip 



uVr, 



A' 



8irs 

x/d(r-r')rf 2 -RVln 



,2 ^dipb) _ s „ 
d p — ^ e 

r-R-f 



R 



V'ln 



T? Ill 



R 



(37) 
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By construction, gradients in Eq. (l37t act on r and 
r'. Equivalently, we can assume that both gradients 
act on R. Integrating by parts and using the identity 
V 2 In |R| = 27T<5(R), we get rid of the R integral. Subse- 
quent integration over r — r' yields 



B 



dip 



2s 



V P A 4 



dm rrC 



z' dip (p) 

Zo 



(38) 



Let us now analyze the determinant factor Z' dip /Z - 
The major part of the excitations near the dipole mini- 
mum of the action have wave lengths in the interval be- 
tween A -1 and m. These modes, living in the relatively 
slow "background field" created by vortices, lead to per- 
turbative renormalization of the parameters entering the 
dipole action S m , in accordance with Eqs. (|2"T1) and (f2"5|) . 
The rest of the determinant provides a non-universal fac- 
tor dependent on the details of the ultraviolet regulariza- 
tion or, equivalently, on the inner structure of the vortex 
core. We do not assume any particular ultraviolet cutoff 
but instead will include all non-universal factors in the 
definition of the vortex fugacity y. Using Eq. (f3"2"j) . we 
represent the dipole contribution to B\ in the following 
form: 



7->dip 

B 1 1 oc 



V p e- 2S *A / dm(Amf~ 2K ^\ 



(39) 



All numerical factors omitted in this expression remain 
finite in the replica limit N — > 0. 

In order to find the dipole contribution to RG equa- 
tions for a and c, we express the derivative of B^ lp in 
terms of fugacity, 



dlnL 



(40) 



All uncontrolled numerical prefactors are now hidden in 
the definition of y: 

y(L) cx ^/Ne^ (AL) 2 ~ K{ * L \ (41) 

Note that the ratio V p /N remains finite in the replica 
limit. Taking derivatives of the Eqs. (|T6| and (|20|) and 
using the identity Q36p , we readily obtain the dipole con- 
tribution to the RG equations in terms of fugacity y. 
Adding the perturbative corrections from Eqs. ([27)1 and 
(f28|) . we get the result 



ay 

dhxL ) ay 2 



l + N/2' 



AIII, BDI, 
CII, 



(42) 



dc 



Nc 



1 + JV' 
oy 2 

2 + N 1 



AIII. BDI, 



CII. 



(43) 




FIG. 1: Schematic RG flow near the metal-insulator transition 
in the symmetry class AIII according to Eqs. (144(1 and (|46[1 . 



Renormalization of fugacity follows from Eq. (|41l) . Dif- 
ferentiating with respect to In L and then taking the limit 
AL — >• 1, we obtain 



dy 



(2-K)y. 



(44) 



This equation is expressed in terms of the stiffness pa- 
rameter K defined in Eq. (|30[) . Renormalization of K 
follows from the other two RG equations (|4"2")l and (|43l) . 
In the replica limit N — > 0, the equations for a and K 
acquire the form 



da 
d\nL 



dK 
dlnL 



— < 



-ay 

( 1 

4 " 

1 o cry 2 

- - 2Ky 2 + 



2Ky 2 , 



1 



16 ' 

2 



-2Ky 2 



ay 



4 ' 



AIII, 

CII, 

BDI. 



(45) 



(46) 



From this result we see that vortices provide negative 
localizing correction to the conductivity. In a good metal- 
lic sample with a 3> 1 and exponentially small (in a) 
fugacity y, vortices are bound in tiny dipoles. The over- 
all negative correction to conductivity remains finite at 
long length scales since the fugacity rapidly decreases in 
the course of renormalization. With lowering the start- 
ing value of conductivity and increasing fugacity, larger 
dipoles come into play and negative correction to a is 
more pronounced. When both a and y are of order 1, 
the phase transition occurs. Our RG equations are not 
quantitatively accurate in this limit but can be used for 
a qualitative description of the transition. 



C. Vortex-driven localization 

Let us study in more detail the RG equations in class 
AIII. As was mentioned in the end of Sec. IIV1 Eqs. (|4"5)) 
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and (|46)l are exact to all orders in the 1 /a loop expansion. 
The only assumption used in deriving these equations is 
y <C 1. Let us neglect this assumption for a moment and 
treat the RG equations as if they were exact. Another 
feature of the symmetry class AIII is the partial separa- 
tion of variables. Renormalization of K and y is inde- 
pendent of er. This allows us to study the 2D flow in the 
K - y plane. The flow is depicted in Fig.[TJ The diagram 
contains one unstable fixed point at the position K = 2, 
y = 1/4. The black line passing through this point sep- 
arates the two stable phases of the system. Below this 
line, the flow is directed towards y — and K — > oo. This 
part of the diagram corresponds to metallic phase. With 
vanishing fugacity, renormalization of conductivity stops 
at some finite value, while the stiffness K keeps increas- 
ing solely due to renormalization of the Gade parameter 
c. In the metallic phase, the system rapidly approaches 
the regime described by the Gade sigma model without 
vortices. The finite contribution of the dipoles can be 
included into bare parameters of this sigma model. 

Above the critical line, fugacity y increases indefinitely, 
which signifies vortex unbinding. This phase is analogous 
to the high-temperature disordered phase in the BKT 
problem and corresponds to the insulating state of the 
sample. Indeed, according to Eq. (jH)) , the conductivity 
flows all the way to zero with growing y. Perturbative 
RG equations rapidly become inaccurate in this limit and 
cannot be used deep in the insulating phase. Neverthe- 
less, using analogy with the BKT transition, we conclude 
that the distinction between metal and insulator states 
is robust. 

RG equations (|44|) . (|45|) . and (j46|) were derived as- 
suming small value of the fugacity (and large value of 
cr in classes CII and BDI). This means that the metal- 
insulator transition fixed point lies at the border of the 
applicability region of these RG equations. The theory 
is not quantitatively controllable in the vicinity of the 
fixed point but the overall qualitative picture of the RG 
flow and, in particular, the very presence of the unstable 
fixed point should be correct. Numerical simulations of 
electron transport in 2D chiral disordered systems would 
provide a quantitative estimate of the critical properties 
of the localization transition. 

We have found a fixed point of the metal-insulator 
transition at a — 0. This signifies the absence of mini- 
mal metallic conductivity in the considered problem. In 
a realistic disordered chiral system, the bare value of fu- 
gacity is related to the Drude conductivity by an expo- 
nential law lny oc —a. This means that the critical state 
is achieved when both y and a are of order unity. Hence 
the Anderson transition will normally happen at a ~ 1. 
The exact critical value of a depends on microscopic de- 
tails of a particular realization of the disordered system. 

In the next section we present an alternative frame- 
work for the derivation of non-perturbative renormaliza- 
tion equations for the chiral class AIII. (We expect that 
a generalization of this approach to other chiral classes 
is possible as well but it requires further work.) 



VI. CLASS AIII AND GENERALIZED 
SINE-GORDON THEORY 

A. Sigma-model and duality transformation 

The group nature of the sigma-model target space for 
the class AIII allows us to develop yet another description 
of vortices, bearing analogy with the sine-Gordon theory 
of the BKT transition^ To construct such a description 
we first work out a dual representation of the AIII sigma 
model. The inclusion of vortex excitations turns out to 
be straightforward within this approach. 

As a starting remark, let us note that for any field 
Q(r) of unitary matrices one can define an associated 
two-component vector field 

h=-iQ- 1 VQ. (47) 

The components of h are Hermitian N x N matrices. 
Conversely, given a vector field h with Hermitian entries 
and the value of the unitary matrix Q at a single point, 
one can reconstruct Q in the whole plane by integrating 
the system of equations 

V fl Q = iQh l ,. (48) 

The solution of Eq. (14"5|) . if exists, is automatically a uni- 
tary matrix field. Examining the cross-derivatives of Q, 
we find that the compatibility condition for the system 
(f4*5|) is given by 

F/iv = VJi u - V v h„ + i[hp, K] = 0. (49) 

Quite remarkably, the integrability constraint (|49]) can 
be viewed as a requirement of vanishing field-strength 
tensor for the non-abelian "gauge field" h. 

We thus come to the conclusion that the matrix field 
Q is in one-to-one correspondence with the vector field h 
subject to the integrability constraint (l4"5|) . In particular, 
the grand partition function of the sigma model can be 
represented as 




where the action functional is quadratic in h, cf . Eq. ([3]) : 
S[h] = |^[aTrh 2 + c(Trh) 2 ]. (51) 

A straightforward inspection shows that the integration 
measure over h in Eq. ([5T)]) is indeed flat. 

The functional <5-function in Eq. (|5T)|) can be resolved 
by the integration over an auxiliary Hermitian matrix 
field 8 

Z = J DhDee- s[h ' @] , (52) 
S[h, 6] = S[h) +i [ d 2 r Tr{QF xy ). (53) 
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Since the action S[h, 9] is quadratic in h, we can easily 
perform the integration over h. As a result, we obtain a 
dual representation of the AIII sigma model in terms of 
the matrix O: 



Z = J Dfi(Q)e 



S[0] 



(54) 



S[Q] 



2ir 



a + Nc 



d 2 r (W ) 



+ —{5 llv + ie» v ) I rf^(V^ Q ).g Qb (e)(V^ 6 ). (55) 
a ' 



[i.e. replacing the constraint — by Eq. (|59|) . This 
yields 



Z = / D^{Q)e- s ^ ' 

^ n=l 



(A 2 yo) r 



x ^ Y[ / d 2 r i d\p i )ey.Y>{2TTir] l {pi\e(Yi)\p l 



Dfi(Q)e 



-s[e]-s v [0] 



(60) 



Here we have introduced a decomposition = 8q/ v N + 

X^^Li -1 S a T a of the matrix over the generators of 
V(N). Further, g ab is a (A 2 - 1) x (A 2 - 1) matrix 
given by 



AT -1 



.9 = (1 - 4ttZ/(7) , Z ab = J2 e cJabc, 



(56) 



where f a t c = —i Tr T a [T b , T c ) are the structure constants 
of the su(N) algebra. The measure of functional integra- 
tion over 9 in (J5^ is Dfi(0) = ]J r det{c?[0(r)]} dO{r). 
Since (g + g T )/2 = gg T can be regarded as the metric 
tensor of the model, we recognize that the integration 
measure is consistent with the metric. 

By construction, the model defined in Eqs. (|54)l . (1551) is 
fully equivalent to the original AIII sigma model. In par- 
ticular, the two theories should obey identically the same 
perturbative renormalization group. In Appendix [A] we 
verify this fact explicitly within a one-loop calculation. 

So far we were completely neglecting the vortex exci- 
tations. To include them into the theory, one must real- 
ize that integrability condition Eq. (14T)1) can be slightly 
relaxed. Indeed, as discussed above, a single vortex cen- 
tered at the origin of the coordinate frame is described 
by the Q matrix 



Q = l + |p)(e*-l)<p|. 



(57) 



Here <f> is the polar angle and \p)(p\ projects onto some 
unit vector \p) in the replica space. The corresponding 
h field has only the azimuthal component = \p)(p\/r 
and the field-strength tensor assumes the form 



F xy (r)=2n5(r)\p)(p\. 



(58) 



A straightforward generalization of Eq. (|58p to the case 
of an arbitrary collection of vortices and anti- vortices lo- 
cated at points Vi and characterized by projectors 
reads 

F xy {v) = 27r5^» 7i *(r-r < )bi>(p < |, m = ±1. (59) 



We now repeat the treatment leading to the construc- 
tion of the dual sigma-model representation (|5"4")) . ([55)1 
but this time taking into account vortex configurations 



Here y = e~ Sv is the bare value of fugacity, integration 
with respect to \pi) is performed over all complex vectors 
of unit length with natural integration measure, and the 
vortex-induced correction to the action reads 



2yK 2 / d 2 rd\p) cos27r(p|e(r)|p). (61) 



S y [Q] 



The action of the dual theory given by a sum of Eqs. (|5"5")) 
and (|6Tj) constitutes the central result of this section. It 
provides a convenient starting point for the generation of 
the RG equations (including non-perturbative contribu- 
tions) for the sigma model of class AIII. 



B. RG analysis 

We are now in a position to explore the RG flow of the 
AIII sigma model from the point of view of the dual rep- 
resentation. In the minimal model of class AIII, the U(l) 
sigma model, Eqs. (|55|) and (jSTj) reduce to the standard 
sine- Gordon action 



S[Q] 



S y {0] = 



d 2 7 



2n 



-(V6» ) 2 - 2yA 2 cos27t6> c 



and we recover the BKT renormalization group. Let us 
now explore the contribution of vortices to RG equations 
in a general situation N ^ 1. As in the BKT theory, 
we will do this perturbatively in the vortex fugacity yo- 
We neglect for a while the non-linear terms in S[&], i.e. 
replace S[0] by 



S o [0] = 2ir / d 2 r 



(Wo) 2 
a + Nc 



E 



(W Q ) ; 



(63) 



Following the standard procedure we decompose the 
fields & into the fast and slow modes 6 = 8< + 9>, 
with the fast modes 0> populating a thin shell in mo- 
mentum space A' = A — A A < \k\ < A. The action 
functional becomes 



s[e> + e<] = s [e<] + 5 [e>] + s„[e<] 



5i[e<,e>], 

(64) 
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where 

Si[e>,e<] = 5 c [e>,e<] + 5,[e>,e<], (65) 

^ c(s )[e>,e<] = -2A 2 2/ J d 2 rd\p)U c[s) {v 1 \p)), (66) 

U c (r,\p))= (coB27r(p|e>(r)|p)-l) 

xcos2^(p|e<(r)|p), (67) 

U s (r,\p)) = -sin2 7 r(p|e>(r)|p) sin27r(p|e<(r)|p). 

(68) 

It is easy to see that to the first order in y the inter- 
action term S\ generates the correction to fugacity itself 
governed by the RG equation [the stiffness parameter K 
was defined in Eq. (|30|) ] : 



dy 



<91nA/A 



7 = (2 - K)y. 



(69) 



We now turn to terms of the second order in fugac- 
ity. We observe that, strictly speaking, the action of our 
model does not preserve its form under the RG transfor- 
mations. For example, let us consider the contribution 
to the action due to the interaction S c : 

A5 C [8<] = -\((S 2 )) = -2AV / d 2 Rd 2 rd\ Pl )d\p 2 ) 

x (([/ c (R+,|p 1 ))C/ c (R_,|p 2 ))», (70) 

where R± = R ± r/2. The gradient expansion of the 
slow fields leads to 

AS C [8<] = -AY / d 2 Rd 2 rd\p 1 )d\ P2 ) (X+ + X-), 

x ((cos2 7 r(p 1 |e > (R+)|pi)cos2^(p 2 |e > (R_)|p 2 ))), 

(71) 

with 



X ± = coB27r[(p 1 |e<(R)|pi> ± (p2|e<(R)|pa)] 
1 ± ^-< P i|ve<(R)|p 1 >{p 2 |Ve < (R)b a ) 



(72) 



We see that the RG flow generates terms of the form 
cos27r[(p 1 |6 < |pi)±(p 2 |0 < |p 2 )] corresponding to the cre- 
ation of two vortices (for plus sign) or a vortex and anti- 
vortex (for minus sign) sitting at the same point and char- 
acterized by the projectors |pi)(pi| and |p 2 )(p 2 |. While 
the former process is irrelevant in RG sense, the latter 
one can be even more relevant perturbation than the 
initial cos27r(p|0<|p) term if vectors \pi) and |p 2 ) are 
sufficiently close. However, these terms are suppressed 
by additional power of fugacity y <C 1 as compared to 
the original vortex term, and we neglect them. Clearly, 
this neglect is not justified in the region y ~ 1 where 
the fixed point governing the transition is located. At 
the same time, the theory we are developing is quantita- 
tively controllable only at y <C 1. Discarding the above 



terms should not lead to any qualitative changes in the 
RG flow. 

We are only interested in the most important y 2 terms, 
i.e., those that produce contributions to rcnormalization 
of a and K and are thus responsible for the localization 
transition. To obtain them, we approximate the correc- 
tion (|7"Tj) to the action of the slow fields by 



A5 C [6<] = 2irnA 4 y 2 J d 2 Rd\p) (p|Ve<( J R)|p) 2 
d 2 rr 2 ((cos27r(p|e>(R + )|p) cos27r(p|e>(R_)|p))). 



(73) 



Here k is a numerical coefficient. Averaging now over the 
fast fields and singling out the contributions of the first 
order in AA we get 

A5 C [8<] = 2irK^K 2 y 2 J d 2 Rd\p) (p|Ve<(R)|p) 2 . 

^ (74) 
(We have absorbed an additional numerical coefficient 
into k.) Noting that for arbitrary matrices A and B 



d\ P )(p\A\p)(p\B\p} = 



we finally get 



(75) 



A5 c [e< 



27TK- 



AA V p 
~K~ N(N+ 1) 



K 2 y 2 



(n + l)(V^) 2 + £(V^) 2 



(76) 



Comparing now Eq. ([76]) with the Gaussian action (|63p . 
we conclude that integration over the fast modes has gen- 
erated corrections to a and c given by 



Act = — k 
Ac = — K 



y 2 K 2 o 2 AA 
r(iV + 2) A : 
y 2 K 2 AA 
r(7V + l) A [N+ 1 



2crc + iVc 2 



(77) 



(78) 



So far we were neglecting non-Gaussian terms in the 
action S[&]. When taken into account, they will induce 
the perturbative renormalization of the parameters of the 
model and, in particular, lead to the appearance of the 
scale-dependent a and c in Eqs. (|77| and (f78j). Apart 
from this, the higher order terms do not affect our anal- 
ysis: a straightforward calculation shows that in our ap- 
proximation the interaction term AS S does not produce 
corrections of order y 2 . Combining Eqs. (J77J), and 

(1751 . we obtain a system of RG equations for the AIII 
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sigma model 
da 



-N ■ 



y 2 K 2 a 2 



dlnL " T(N + 2)' 

8K 1 - N y 2 K 2 
dhiL 

dy 



(79) 



T{N + 1) \N + 1 



2a 2 



+ 2ac + Nc' 



dlnL 



(2 - K)y. 



(80) 
(81) 



These equations are equivalent to those obtained in Sec. 
|V]up to subleading terms. To make this equivalence ap- 
parent, we redefine the fugacity parameter according to 
V ^ (y/K)y / T(N + \)ja. This amounts to changing 
the (uncontrolled) pre-exponential factor in the exponen- 
tially small quantity y. After such a rescaling, the RG 
equations (|42l) . (l43l) . and (l44l) are reproduced. 



VII. TOPOLOGICAL TERMS 

In previous sections, we discussed localization within 
the Gade-Wegner sigma model, Eq. ([3]), of three chi- 
ral classes. Now we will consider the situations when 
the sigma-model action is augmented with an additional 
topological term. The topology of the sigma-model tar- 
get spaces allows (in two dimensions) inclusion of the 
Wess-Zumino term in class AIII and the Z 2 term in 
class CII. We will demonstrate below that these extra 
terms (arising in models of random Dirac fermions of 
the corresponding symmetries) crucially affect localiza- 
tion properties. Specifically, we show that vortex excita- 
tions do not appear when a topological term is present. 



A. Class AIII with Wess-Zumino term 

The field of the sigma model of class AIII is a uni- 
tary matrix Q(r) € \J(N), with r belonging to the 2D 
coordinate space. In our analysis we consider only field 
configurations with the Q matrix taking some fixed value 
at spatial infinity. (Otherwise the sigma-model action in- 
evitably diverges due to gradient terms.) This allows us 
to compactify the coordinate space making it equivalent 
to a two-sphere. The compactified real space S 2 can be 
viewed as a surface of a three-dimensional solid ball and 
we introduce the radial coordinate r such that r = 1 at 
the surface (i.e., in the physical 2D space) and r = in 
the center of the ball. The matrix Q can be continuously 
extended to the interior of the ball such that 



Q(r,r) 




(82) 



This extension is always possible since the second homo- 
topy group of the target space is trivial, 7r 2 (C/) = 0. 



In terms of the extended matrix Q(t), the Wess- 
Zumino term acquires the form 



& 



ik 



WZlV] - T7T- ^ 
127T 

x J dT(PrTr(Q- 1 V li QQ- 1 V u QQ- 1 V x Q). (83) 

The integrand in this expression explicitly depends on 
the values of Q at r 7^ 1, i.e., away from the physical 2D 
space. However, the variation of the Wess-Zumino action 
can be represented as an integral of a three-dimensional 
vector divergence: 



8Sy/z[Q] = -7—£uv\ 

47T 



X 

ik 
4tt 



J dr rf 2 rV M Tr (q^SQQ-^uQQ-^xQ) 



'Tr(Q- 1 8Q[Q- 1 V x Q,Q- 1 V v Q\y (84) 



Thus the actual value of the Wess-Zumino term is deter- 
mined only by the physical values of Q at r = 1 up to a 
constant. This constant does not change with small vari- 
ations of Q away from the physical 2D space but takes 
different values for topologically distinct extensions of Q 
in the third dimension. These nonequivalent extensions 
are classified by the third homotopy group ^(t/) = Z. 
For any two extensions the values of the Wess-Zumino 
term differ by 2itik times an integer number. Thus the 
Wess-Zumino theory is well-defined for any integer value 
of A. 

Introducing vortices in the Wess-Zumino theory is 
problematic. In order to avoid the singularity in the cen- 
ter of a vortex, we exclude a small region of the vortex 
core from our physical space. This introduces a bound- 
ary in the problem and Q is not a constant along this 
boundary. As a result, the 2D physical space cannot be 
compactified to the two-sphere. Thus the construction of 
the Wess-Zumino term, involving extension to the third 
dimension, becomes ill-defined. 

One naive way to overcome this difficulty is to use a 
local 2D representation of the Wess-Zumino term4£ As 
was discussed above, the Wess-Zumino term actually de- 
pends only on the values of Q in the physical space. Us- 
ing any explicit parameterization of the unitary matrix 
Q by a set of coordinates ipi, the Wess-Zumino term can 
be written as a 2D integral of a suitable skew-symmetric 
differential form A,-,-: 



S- 



wz 



= ik 1 



d 2 r A* 



(r)] V^V^j 



(85) 



In this local representation, one can integrate the differ- 
ential form over a finite physical space with boundary. 
Such a construction is however unsatisfactory because 
it violates the global gauge symmetry of the system. In- 
deed, the Wess-Zumino action (183)) is manifestly invariant 
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under the global transformation Q i— > UlQUr parame- 
terized by two constant unitary matrices Ul,r- However, 
the local density of the Wess-Zumino term, Eq. (|85|) . is 
not invariant under global rotation of fields. Instead, the 



tensor Xij transforms as 



X ij n- X ij 



dtpi 



(86) 



where a set of functions f3(ip) encodes the information 
about Ul and Ur. The local expression for the Wess- 
Zumino term changes by an integral of a total derivative: 



Swz Swz + 2ike flL , J d 2 r V Al (/3 J -V w ipj). 



(87) 



If the model is considered on a manifold without a bound- 
ary, the above integral vanishes and the Wess-Zumino 
term is indeed invariant under the gauge transformation. 
If, however, the real space integration is performed over 
a bounded region, the integral in Eq. (j5Tf yields the cir- 
culation of f3j\Jipj along the boundary and may become 
non-zero. This signifies the breakdown of the gauge sym- 
metry at the boundary. 

In order to understand better the boundary effects in 
the Wess-Zumino theory, we will resort to the disordered 
fermion problem yielding the class AIII sigma model with 
the Wess-Zumino term. The typical example is given by 
disordered massless Dirac fermions with random vector 
potential. The Hamiltonian of such a model has the form 



H 



.Pa 





i P y 



Px ~iPy+' 



Here a is a random complex- valued matrix acting in the 
auxiliary flavor space. The only symmetry of the Hamil- 
tonian is the chiral symmetry H = —a z Ha z . Since the 
spectrum of the Hamiltonian is unbounded, it appears 
to be impossible to introduce the boundary condition (a 
classically forbidden region for massless Dirac fermions) 
preserving the chiral symmetry. In order to model a 
hole in the sample, one has to add an additional mass 
term ma z to the Hamiltonian and consider the limit of 
large m. Such a term explicitly breaks the chiral sym- 
metry. In the corresponding sigma model, the Q matrix 
at the hole boundary will be restricted to the manifold 
M a = U(N) /V(N/2) x XJ(N/2) of class A. Such a bound- 
ary condition will maintain only the diagonal part of the 
global gauge symmetry Ul = . This also forbids 
the vortex excitations inside the hole since Ma is simply 
connected. Thus we see that the short distance regu- 
larization (making a small hole), needed to introduce a 
vortex in the sigma model, breaks the chiral symmetry 
and does not allow a vortex excitation. We conclude that 
vortices are incompatible with the Wess-Zumino term in 
the action of class AIII. 

The model of massless Dirac fermions in a random 
magnetic field, described by the Hamiltonian (|88| is ex- 
actly solvable. The coupling constant characterizing the 



vector potential strength is exactly marginal, so that the 
model possess a line of fixed points. A remarkable prop- 
erty of this problem is that the system never gets truly lo- 
calized, however strong the disorder is. In particular, the 
conductivity is equal to e 2 /irh for any disorder. There- 
fore, the absence of the localizing vortex contribution 
in the corresponding sigma-model is consistent with the 
known exact solution of the underlying fermionic prob- 
lem. 

Let us illustrate the incompatibility of vortices and 
Wess-Zumino term in a more explicit way. We will con- 
struct an extension of the sigma-model manifold such 
that the theory will be well-defined inside the vortex core 
while all the symmetries are preserved. Let us extend the 
Q matrix by one row and one column embedding U(iV) 
into SU(A + 1). We will associate a mass M with the 
extra degrees of freedom suppressing them away from a 
vortex core. Explicitly, consider the following action for 
the extended unitary matrix Q: 



S[Q] 



87 



(jTriQ^WQ) 2 + M 2 TriQRQ- 1 R) 



Here the matrix R of the form 



R 




Swz[Q]- (89) 



(90) 




is introduced to single out the off-diagonal elements in 
the first row and first column of Q. These elements are 
made massive by the second term of the action (|89l) . 

Away from vortices the matrix Q contains only the soft 
modes arranged as follows: 



(91) 



The lower right diagonal block is nothing but the unitary 
matrix Q while the upper left diagonal element is fixed 
such that detQ = 1. With such a form of Q, the action 
(f89|) coincides with the standard AIII class sigma-model 
action with a Wess-Zumino term added. At the same 
time, the group SU(iV + 1) is simply connected hence 
the vortices are topologically trivial configurations in the 
extended model. Inside a vortex core, massive elements 
of Q become non-zero and provide an overall smooth field 
configuration. The size of the core is determined by the 
competition between energy loss due to the mass M and 
energy gain due to avoiding large gradients. This yields 
the core size ~ \fa/M. 

Within the extended model we can examine the inner 
structure of the vortex core. Assume for simplicity that 
the vortex occurs in the first replica, i.e., the correspond- 
ing unit vector is \p) = {1, 0, . . . , 0}. The whole vortex 
configuration, including the core, will involve only the 
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upper- left 2x2 block of the matrix Q. This block is an 
SU(2) matrix and we can explicitly parameterize it by 
three angles in the following way: 



ism9e lx cosOe 1 ^ 
\ 



1 






JV-lxJV-1, 



(92) 



Far from the vortex core, the angle 9 vanishes and the 
matrix Q acquires the form (|9ip ; it is independent of x- 
Going around the vortex, angle (f> rotates by 27r. In the 
center of the vortex 9 = Tt/2 and Q is independent of <f> 
but explicitly depends on \- Using the ansatz (|92p . we 
can minimize the action (|89|) . The symmetry of vortex 
allows us to fix parameter </> equal to the polar angle 
and x constant. We see that the vortex core acquires 
an inner U(l) degree of freedom, x, which is beyond the 
sigma model and is effective only in the extended theory. 
The action is minimized by a proper 9(r) dependence. 

The Wess-Zumino term is responsible for the imagi- 
nary part of the action. We can calculate this imaginary 
part without solving for 9{r) dependence. Let us consider 
the variation of the Wess-Zumino term with respect to 
spatially constant x- Upon substitution of Eq. (|92[) into 
Eq. ([SI]), we obtain 

SS WZ [Q] = ^ J d 2 r sm29(V x 9V y 4>~V y 9V x 4>) 
= ikSx sin 2 9(r = oo) - sin 2 9(r = 0) = -ikSx- (93) 

Thus the value of the Wess-Zumino term explicitly de- 
pends on x- 

The imaginary part of the action makes vortex fugac- 
ity complex. The arameter x-, that determines the phase 
of the complex fugacity, is an internal degree of freedom 
of the vortex core. Calculating the partition function of 
the system, we have to integrate over x f° r each vor- 
tex. Once the action contains the Wess-Zumino term, 
i.e., k 0, such an integration will exactly cancel the 
statistical weight of each vortex making fugacity effec- 
tively zero. This once again demonstrates that the Wess- 
Zumino model of class AIII does not allow vortex exci- 
tations. 



B. Class CII with Z 2 term 

Let us now consider the system of symmetry class 
CII. The topology of the target manifold admits the Z2 
topological term. It is related to the homotopy group 
iT2(U(N)/0(N)) = Z2. One particular realization of the 
symmetry class CII is provided by the disordered mass- 
less Dirac Hamiltonian of the form Eq. ([T]) with 



Px - ip y 



Px + iPy. 



(94) 



Here a is a random n x n matrix with complex entries. 
The block h obeys symplectic constraint h* = T y hT y . In 
other words, h is an n x n matrix of real quaternions. 
Thus the Hamiltonian built from h indeed belongs to the 
symmetry class CII. 

The derivation of the sigma model for the disordered 
system described by the Hamiltonian ([1]), (|94| is outlined 
in Appcndix[Bj The action of the model has the standard 
form ([3]) with an additional Z2 topological term when n 
is odd. This topological term can be expressed in the 
form very similar to the Wess-Zumino term (|83[) . Specif- 
ically, we have to continuously extend the Q matrix to 
the auxiliary third dimension r according to Eq. (|82[) . 
Then the topological term can be written in the form of 
Eq. (1831) with k = n. Since the second homotopy group 
of the target manifold is non-trivial, the extension (|82|) 
is not always possible within the target space of class 
CII. In fact, it is only possible for topologically trivial 
configurations of Q(r). To apply the Wess-Zumino con- 
struction to a general field configuration, we will assume 
that away from the physical 2D space, at r 7^ 1, Q is any 
unrestricted unitary matrix from U(-ZV), while for r = 1 
it is unitary and symmetric hence belongs to the coset 
space V(N)/0(N) of the class CII. 

Similarly to the class AIII, the value of the Wess- 
Zumino term is actually determined by the physical part 
of Q at r = 1 which is a unitary symmetric matrix. The 
variation of the Wess-Zumino term with small variations 
of Q is given by Eq. (f84|) . Using the property Q = Q T 
and transposing the argument of the trace in the last line 
of Eq. (|M)) . we see that the variation changes sign and 
hence is identically zero. This proves that such a Wess- 
Zumino term in the class CII possesses the main property 
of the theta term: it depends only on the topology of the 
field configuration. 

Apart from the class of topologically trivial configura- 
tions, there is only one extra non-trivial class. In other 
words, the homotopy group ~ki = TLi implies existence of 
localized topological "excitations" , %2 instantons. They 
are their own "antiparticles" ; configuration of two such 
instantons is topologically trivial, i.e. the instantons can 
be brought close to each other and annihilated by an ap- 
propriate continuous transformation of Q. In order to 
prove that the Wess-Zumino term (|53")) . being constant 
in each topological class, distinguishes between them, it 
suffices to evaluate it for one particular non-trivial in- 
stanton configuration. 

Let us consider the minimal model Q £ U(2)/0(2). 
(In fact, in this case the homotopy group 7T2 = Z is 
readier than in the general case N > 2 and the the- 
ory possesses usual Z instantons similar to, e.g., 0(3) 
vector sigma model.) In fact, it is sufficient to consider 
an even smaller target space SU(2)/0(2) since the deter- 
minant of Q anyway drops from the Wess-Zumino term 
We parametrize the SU(2) matrix by three angles 
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in the following way: 



Q = 



cos a cos \e 



i sin a cos \ + sin \ 



i sm v cos x — sin x cost) cos %e 



(95) 



The symmetry condition Q = Q T hxes x = in the 
physical 2D space. For the instanton, we can assume 
that 4> is equal to the polar angle while 9 depends only 
on the radial coordinate and changes continuously from 
— 7r/2 in the center of the instanton to 7r/2 at infinity. 
Extending to the third dimension, we will assume x to 
change from at r = 1 to either n/2 or — n/2 at r = 0. 
Any of these two values uniquely fixes the whole matrix 
Q(r = 0). 

Wess-Zumino action (|83p can be explicitly written in 
the parametrization (|95[) as 

<Swz[Q] = — e^A / d 2 rdr cos 6» cos 2 xV^^V^c/) V A X- 

(96) 

For the instanton configuration, and x depend on r and 
r, respectively, while is the polar angle. Calculating the 
integral, we obtain 



'wz 



ink r 



oo 



i0(r = 0) 



^fink. 
(97) 

The two signs in this expression correspond to the two 
extensions x( T = 0) = ±7r/2. For an odd value of k, the 
instanton action acquires a non-trivial imaginary contri- 
bution. Thus the Wess-Zumino term indeed plays the 
role of a theta term yielding ink times the topological 
charge of the field configuration. 

Once the explicit form of the Z2 topological term is es- 
tablished, we can discuss its interplay with vortices. We 
have already argued that the Wess-Zumino term in the 
sigma model of class AIII makes vortex excitations in- 
effective. Similar arguments can be applied to the class 
CII with the topological term since it has the same struc- 
ture as the Wess-Zumino term. The only difference in 
the class CII is an additional constraint related to the 
time-reversal symmetry. Namely, the statistical weight 
of any field configuration of the class CII sigma model 
must be real. Equivalently, the imaginary part of the 
sigma-model action must be an integer multiple of in. It 
is the topological term that provides this imaginary part. 
Consider an extension of the model from \](N)/0(N) up 
to SU(A + 1)/0(N + 1). Such an extension is given by, 
e.g., Eq. (|89l) with an additional symmetry constraint 
Q = Q T . The vortex configuration in this extended 
model has the form Eq. with the angle x taking 
either or n value. Thus the internal U(l) parameter 
associated with the vortex in the symmetry class AIII 
becomes a I2 degree of freedom in the class CII. The 
two vortex configurations with x = and x — 77 differ by 
an in term in the action, cf. Eq. (1931) . Thus summation 
over the two values of this internal degree of freedom ef- 
fectively annihilates vortex contribution to the partition 
function of the system. 



The above discussion is based on the specific form of 
the extended model (|59")l describing the vortex core. In 
fact, we can lift this restriction and show that a vortex 
possesses an internal Z2 degree of freedom without as- 
suming any particular structure of its core. Consider a 
vortex configuration in the class CII sigma model. We 
suppose that some additional massive degrees of freedom 
become relevant in the center of the vortex and once they 
are taken into account the field configuration is continu- 
ous. The extended model must possess the time-reversal 
symmetry characteristic for the class CII. Hence the sta- 
tistical weight in the extended model is real and the ac- 
tion is real up to an integer multiple of in. Let us now 
create a small instanton far from the vortex center. This 
will change the imaginary part of the action by n. Within 
the extended theory, we can bring the instanton close to 
the vortex center and "hide" it inside the core by a con- 
tinuous field transformation. Since the imaginary part 
of the action takes only discrete values, this transforma- 
tion will not remove an extra in from the action related 
to the instanton. We have thus demonstrated the ex- 
istence of two topologically distinct vortex solutions in 
the extended theory with opposite signs of their statis- 
tical weights. Since there is no other general distinction 
between these two solutions, the real parts of their ac- 
tion must be equal. This once again shows that the total 
statistical weight of a vortex configuration is zero if the 
underlying sigma model contains the Z2 topological term. 

To conclude, additional topological terms in the sigma 
model of both AIII and CII symmetry classes suppress 
formation of vortices and thus prevent the system from 
localization. 



VIII. SUMMARY AND OUTLOOK 

In this paper, we have developed a field-theoretical 
(sigma-model) approach to Anderson localization in 2D 
disordered systems of chiral symmetry classes (AIII, BDI, 
CII). A remarkable feature of sigma models for these 
classes is that the quantum interference effects leading 
to renormalization of conductivity (and thus to Ander- 
son localization) are absent to all orders of perturba- 
tion theory. We have shown that Anderson localiza- 
tion does exist within these models and is governed by 
a non-perturbative mechanism. Specifically, the local- 
ization is due to topological excitations - vortices - of 
the sigma model field. We have derived the correspond- 
ing renormalization group equations which include non- 
perturbative contributions. Analyzing them, we find 
that the 2D disordered systems of chiral classes undergo 
a metal-insulator transition driven by topologically in- 
duced Anderson localization. 

While the mechanism of the localization transition 
- proliferation of vortices — bears an analogy with 
the Berezinskii-Kosterlitz-Thouless transition in systems 
with U(l) symmetry, our RG equations are essentially 
different. The reason for this is a more complex struc- 
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ture of the theory: it is characterized by three coupling 
constants (conductivity a, Gade coupling c, and fugacity 
y) instead of two couplings of the BKT transition theory 
(spin stiffness and fugacity). As a result, the fixed point 
governing the transition turns out to be at non-zero fu- 
gacity. For this reason, the critical behavior at the tran- 
sition cannot be determined in a controllable way. The 
one-loop analysis suggests that this behavior is of power- 
law type (i.e. is more similar to the critical behavior at 
Anderson transitions in conventional classes rather than 
to that at BKT transition). 

For the chiral unitary class AIII, we have presented an 
alternative derivation of the renormalization group based 
on a mapping of the sigma model onto a dual theory. The 
latter has the form of a generalized sine-Gordon theory. 

We have also considered 2D disordered systems formed 
on surfaces of 3D topological insulators of chiral symme- 
try classes AIII and GIL In this case the sigma model is 
supplemented by a term of topological origin: the Wess- 
Zumino term for the class AIII and Z2 theta term for the 
class CII. We have shown that such terms overpower the 
effect of vortices, thus ensuring the protection of surface 
states against the vortex-induced Anderson localization. 

Our work opens perspectives for research in a number 
of important directions. Below, we briefly discuss several 
of them. 

First, it would be very interesting to investigate the 
metal-insulator transition in chiral classes and the as- 
sociated critical behavior numerically. Remarkably, this 
issue is almost unexplored by now. This is in stark con- 
trast with conventional symmetry classes (i.e., symplectic 
class metal-insulator transition and unitary class quan- 
tum Hall transition in 2D, orthogonal class transition in 
3D etc.) where very detailed studies have been carried 
out. Since we are dealing here with non-interacting sys- 
tems in a relatively low (2D) dimensionality, a sufficiently 
accurate numerical analysis is expected to be feasible. It 
would be also interesting to simulate the sigma model di- 
rectly. This would allow to verify the importance of topo- 
logical excitations for localization and to test our pre- 
dictions. Such an approach can be implemented within 
supersymmetric version of the sigma model with Grass- 
mann degrees of freedom integrated out4^ 

Second, it remains to develop a dual, sine-Gordon-like 
theory of the transition, analogous to that presented in 
Sec. IVII for class AIII, for the other two chiral classes. 
Furthermore, we feel that geometric aspects of sigma- 
model renormalization within this dual formalism deserve 
a more thorough investigation. 

Third, a natural question arises concerning the metal- 
insulator transitions in chiral classes in 3D (and higher 
dimensionalities). We expect that also there the tran- 
sition will be driven by topological excitations, namely, 
vortex lines. 

Fourth, in analogy with Z vortices studied above, 2D 
sigma models of two classes, All and Dili, allow for Z2 
vortices. The difference is that these two classes do show 
quantum interference effects on the perturbative level. 



Therefore, in contrast to the chiral classes, the vortices 
in classes All and Dili will not constitute the only driv- 
ing mechanism of Anderson localization but rather will 
contribute to renormalization of conductivity along with 
perturbative terms. 



Fifth, interaction effects play an important role in low- 
dimensional systems and may strongly affect the nature 
of the metal-insulator transition. Two-dimensional Dirac 
fermions subjected to short-range interaction exhibit the 
Mott transition 42 unlike the non-interacting case when 
the system remains metallic. It would be very interesting 
to investigate the interplay of interaction and vortices in 
systems with chiral symmetry. 



Finally, we close with a more general comment. The 
importance of topological aspects of field theories of dis- 
ordered systems was recently emphasized in the con- 
text of topological insulators 16,21 . There, a possibility of 
emergence of a Wess-Zumino term or Z2 theta term in the 
corresponding dimensionality and symmetry class signals 
the existence of a topological insulator phase. There is 
at present a growing appreciation of the fact that the 
topological properties of sigma- model manifolds may cru- 
cially affect the physical observables even for those com- 
binations of dimensionalities and symmetries than do not 
allow for topological insulators. In particular, a recent 
work 43 has shown that a particular symmetry of local 
density of states distributions holds for sigma models 
(and thus for critical points) of five symmetry classes (A, 
AI, All, C, CI) but dos not hold for the remaining five 
(AIII, BDI, CII, D, Dili). The distinct feature of the lat- 
ter five classes is the presence of U(l) or 0(1) subgroup in 
the sigma model target space M leading to a non-trivial 
topology: 7Ti(A4) = Z for chiral classes (AIII, BDI, CII) 
and 7To(A4) — Z2 for Bogoliubov-de Gennes classes D and 
DHL It is worth emphasizing that these topologies ren- 
der these five classes topological insulators in ID. Equiv- 
alently, models of these symmetry classes may support 
eigenstates with exactly zero energies^. As the paper 
Ref. IH showed, the same U(l) and 0(1) degrees of free- 
dom that are responsible for topological insulator prop- 
erties in ID in fact crucially affect the multifractal spec- 
tra at higher dimensionalities. The present work shows 
that the U(l) topology of sigma-model manifolds of chi- 
ral classes is also responsible for Anderson localization in 
these classes in 2D (and likely also in higher dimensions) . 
Earlier, the authors of Ref. \44 argued that the 0(1) = Z2 
degree of freedom of the sigma model is responsible for 
localization in class D in two dimensions. Thus, it turns 
out that the importance of topological aspects of field 
theories of disordered systems goes well beyond that ex- 
pected on the basis of classification of topological insu- 
lators. Full ramifications of these observations remain to 
be understood. 
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The renormalization of the theory can now be carried 
out in the standard manner. We split the fields 8 a into 
the fast and slow components, 8 = 9 < + 6*>. We should 
then expand in the fast fields up to the second order 
and integrate them out. It proves convenient to perform 
the change of fast variables 6 > = #(#<)#> so that the 
decomposition of 8 into the fast and slow fields reads 



1-Z K 



(A3) 



Appendix A: Perturbative renormalization in dual 
representation 

The purpose of this appendix is to illustrate the equiv- 
alence of the model (|5"4")l . ([55]) to the original AIII 
sigma model by a one-loop perturbative renormalization. 
Within this calculation we can ignore the Gaussian field 
8q in Eq. (|55|) (since renormalization of this sector is 
trivial on the perturbative level) and concentrate on the 
renormalization of the non-abelian sector of the theory 
parametrized by fields 8 a with a — 1 , . . . N 2 — 1 . To 
make more explicit the possibility of a loop expansion 
controlled by the parameter a 3> 1, it is also convenient 
to rescale 8 a by a factor a /Air. The action of the model 
acquires now the form 

Sffl = J^( S M» + ie ^) J d 2 r(y^8 a )g ab {8){V v 8 b ), 

(Al) 

N 2 -l 

g(0) = [1 - Z(9)]-\ Z ab (8) = Y, e cfabc (A2) 

c=l 



Here we introduced for notational brevity Z < = Z{8 < ). 
This change of integration variables cancels the contri- 
bution of the non-trivial integration measure Dfi(9) to 
the renormalization group equations and guarantees the 
absence of linearly diverging diagrams in one-loop calcu- 
lation. Denoting also Z(9>) by Z > we have 



V0=(1 + Z>)V0< + (1-Z<)V0>, (A4) 
Z{8) = Z K + Z> + [Z> , Z<] , (A5) 



This implies that the quadratic-in-#> contribution to the 
action reads Si — + SS2 , 



c(0) _ °" 
b2 "8^ 



d 2 rV6^V6>>, 



SS 2 



l€ L , 



8n 



d 2 r [V P 0£(1 + Z<) + V M 0£<7<Z>Z<] [V„0> + <?<(2 - Z < )Z >g< V ,8^ - S, 



(0) 
2 • 



(A6) 
(A7) 



Performing now the expansion in SS2, averaging over 
fast fluctuations, and retaining only the logarithmically 
diverging contributions, we find the correction to the ac- 
tion functional of the slow fields: 

N A 

x J cf 2 rV^[2# + 30 2 ]V„0< . (A8) 
We can recast the action for the slow fields into the 



original form by correcting the conductivity a <— > a and 
switching to the rescaled filed 

*-'( 1 -7 b1 ^)' (A9) 

»-.(l + £h4). (AlO, 

(Note that, after the rescaling, the factor 1/(1 — Z) should 
be expanded to the first order in 1/tr, which is the ac- 
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curacy of the one-loop calculation.) We see now that 
the dual model defined by Eqs. (jM)) , (|55|) reproduces the 
correct renormalization of the conductivity a. 



Appendix B: Chiral sigma model for massless Dirac 
fermions 



In this Appendix we outline the derivation of the sigma 
models for the massless Dirac Hamiltonians of AIII and 
CII symmetries. Apart from demonstrating the conse- 
quences of the chiral symmetry in the sigma-model lan- 
guage, we will also discuss possible topological terms. 
Such terms frequently appear due to chiral anomaly of 
the massless Dirac electrons. 



1. Class AIII 

Let us first consider the AIII symmetry class. A gen- 
eral chiral Hamiltonian has the block off-diagonal struc- 
ture (UJ . We consider the following Dirac Hamiltonian of 
AIII symmetry [cf. Eq. 



H 




(Bl) 



Here p± = p x ± ip y and a is a random complex- valued 
matrix of size n x n. Further, we adopt the most stan- 
dard model of Gaussian white noise disorder with the 
correlator 



(a*,a k i) = — 5ik5 
J n 



ikOjl- 



(B2) 



Parameter a quantifies the disorder strength and de- 
termines the scattering rate of electrons. Within self- 
consistent Born approximation, the scattering rate at the 
Dirac point (zero chemical potential) is given by 



7 = A e 



-2ir/c 



(B3) 



where A is the effective band width (maximal allowed 
energy) for Dirac Hamiltonian. The self-consistent Born 
approximation is valid in the limit n 3> 1. This also 
corresponds to high Drude conductivity n/ir(e 2 /h) at the 
Dirac point and justifies the applicability of the sigma 
model. 

Derivation of the non-linear sigma model starts with 
the replicated action written in terms of fermionic (anti- 
commuting) fields: 



S : 



o 



a 




(B4) 



The lower indices i and j take n possible values in the 
flavor space while the upper index a enumerate N repli- 
cas. We proceed with averaging e~ s over the random 



matrix a. With the help of Eq. (|B2I) . the effective action 
acquires the form 



S = $fp_W + tfp+tf + 5 dis , 



5,i 



www- 



(B5) 
(B6) 



The quartic term Sdis is further decoupled with the 
help of an auxiliary complex matrix field Q acting in 
replica space. This is achieved by adding to the action 
a term Tr Q^Q and shifting the Q variable by a suitable 
quadratic expression in fermions. 



<3dis h ~?' iJdis n 

a 



u + — c# 

rvy 



(B7) 



— TrQtQ - ilUfQabtf + $Q* ab W 
a V 



Now we add the rest of the action from Eq. (|B5|) and 
then perform Gaussian integration over fermion fields. 



5=^Tr, 



^ I^TrQtQ-nTrln 
a 



-i-yQ 


V 


P+ 




-i'yQ 


v N 




-*7<3 1 ; 



(B8) 



The next step of the sigma model derivation involves 
saddle-point analysis of the above action for Q. The 
saddle-point equation is equivalent to the equation of self- 
consistent Born approximation with ijQ playing the role 
of self energy. Therefore one particular solution is just 
the unit matrix Q = 1. Other solutions can be found 
by unitary rotating the fermion fields in the first line of 
Eq. (|B8|) such that the p± terms remain intact. These 
rotations generate the gauge group U(iV) x U(iV) rotat- 
ing Q by the two unitary matrices from left and right. 
Thus the matrix Q takes values from the symmetric sub- 
group U(N) x U(A)/U(A) = U(N) of the global gauge 
group. This establishes the manifold of the class AIII 
sigma model. 

The sigma-model action is a result of the gradient ex- 
pansion of Eq. (|B8|) with a slowly varying unitary matrix 
Q. This gradient expansion should be carried out with 
care in view of the chiral anomaly of the Dirac oper- 
ator under logarithm. A systematic description of the 
expansion procedure, including the methods to treat the 
anomaly, can be found, in particular, in Ref. Il5l . The 
result of the gradient expansion reads 



S[Q] 



n 
8^ 



Tr(Q+VQ) 2 + S 



wz| 



(B9) 



This is the standard sigma-model action of class AIII, 
Eq. ([3]) with a — n and c — 0, with an additional Wess- 
Zumino term of the level k = n, see Eq. Appearance 
of the Wess-Zumino term is the direct consequence of the 
chiral anomaly of the Dirac Hamiltonian. The structure 
and properties of the Wess-Zumino term are discussed in 
the main text, Sec. IVII Al 
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Class CII 



Now we discuss the derivation of the sigma model with 
Z 2 topological term in the symmetry class CII. Consider 
the Hamiltonian [cf. Eq. ((Ml) ] 



H = 




h = 




(BIO) 



As in the previous case, a is a random complex-valued 
matrix of size n. The block h fulfills the symmetry condi- 
tion h = Tyh*T y and hence represents a (operator- valued) 
real quaternion matrix. Thus the Hamiltonian (|B10j) in- 
deed belongs to the symmetry class CII. 

We assume the same Gaussian white-noise distribu- 
tion of a as in the previous section, Eq. (|B2|) . Disorder- 
induced scattering rate, Eq. ()B3|) . is also reproduced 
within the self-consistent Born approximation; Drude 
conductivity is twice larger, 2n/7r(e 2 //i), since the Hamil- 
tonian (|B10[) contains 2n coupled Dirac fermions. 

As in the previous section, we start with the replicated 
fermionic action, 



= + + - ■!>:' (bh) 

Here we have introduced the doubled fermionic fields 




* = 



$ = 



(B12) 

Starting from the second line of Eq. (|B11[) and further on, 
we assume that the replica index a takes on 2N values 
running through N replicas and two components of the 

vectors (lBl"2|) . 

We average e~ s from Eq. (jBlip using the correlator 
Eq. (IB"2l and obtain 



(B13) 
(B14) 



S = <&?p_*? + $?p+$- + S d 

S dis = - tfJ^S*?. 
n J J 



These expression are very similar to Eqs. (|B5[) and (|B6I) . 
The only difference is in the order of replica indices in 



the quartic term Eq. (|B14|) . As before, we introduce the 
matrix field Q and decouple the action according to 



"Sdis l— ^ ^dis ' 



717 



ia 



Qab + —* b M 



Qab 



ia 



717 



TV) 



Tr Q^Q-i^Qab^ + ^Ql^t)- (B15) 



Adding the clean part of the action from Eq. (|B13I) and 
performing Gaussian integration over fermion fields, we 
obtain the action in the form 



h-> — TrQtQ-nTrln 
a 



t-ijQ 


P- 


V p+ 




-ilQ 


p - ) 


. P+ - 


-ilQ*) 



(B16) 



This result is again very similar with the class AIII ex- 
pression (|B8|) . The only difference is that the matrix Q' 
under the logarithm in Eq. (|B8I) is replaced with Q* in 
the present case. The standard saddle point of the above 
action is Q = 1. Other saddle points are generated by ro- 
tations of the fermion fields in the first line of Eq. (|B16[) . 
As follows from the saddle point analysis, these rota- 
tions should maintain the complex conjugacy of Q and 
Q* in the argument of the logarithm. This is achieved 
by the unitary transformation of the form $ t—¥ J7$, 
* 1-^ U*^>, $ 3>C/ t , § h> vj/[/ T parametrized by the 
unitary matrix U of the size 2N. The global gauge group 
is thus U(2N). The saddle manifold is parametrized as 
Q = U T U and contains all symmetric unitary matrices, 
Q e U(2A)/0(2A), as it should be for the sigma model 
of class CII. 



For any symmetric unitary matrix Q, the action (|B16p 
is identical to Eq. (|B8[) . This means that the sigma- model 
action Eq. (|B9[) obtained by the gradient expansion of Eq. 
(|B8|) , is also valid for the symmetry class CII provided the 
constraint Q = Q T is maintained. With this restriction, 
the Wess-Zumino term of the level k = n becomes a Z2 
topological term for any odd n, as we discuss in Sec. 
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